Radiation can never again dominate Matter in a Vacuum Dominated Universe 
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We demonstrate that in a vacuum-energy-dominated expansion phase, surprisingly neither the 
decay of matter nor matter-antimatter annihilation into relativistic particles can ever cause radiation 
to once again dominate over matter in the future history of the universe. 
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The study of decaying particles in a cosmological con- 
text has a long history % 0, 0, H i, 0, 0, I, ©, BjJ. 
Although there are many variations on this theme, there 
is one constant result: if the universe is initially domi- 
nated by nonrelativistic particles, and this matter under- 
goes a standard exponential decay into relativistic parti- 
cles, then the universe rapidly transitions from a matter- 
dominated phase into a radiation-dominated phase when 
the age of the universe is roughly equal to the particle 
lifetime. 

In this paper we show, rather suprisingly, that this con- 
clusion is no longer valid at the present time, if the dark 
energy reflects a constant vacuum energy density in the 
universe. Once the universe begins to enter a phase of 
exponential expansion, no process leading to the disap- 
pearance of matter can ever cause radiation to once again 
dominate over the remnant matter density, no matter 
how small. 

Consider a nonrelativistic component with density pm, 
decaying into a relativistic component pn, and take the 
lifetime of this decay to be r, with a standard exponen- 
tial decay law. (For simplicity, we will assume that the 
decay products are "sterile", i.e., that they do not inter- 
act significantly with anything else, but our conclusions 
do not depend on this assumption). Then the equations 
governing the matter and radiation evolution are 

dpM „ u I 

= SHpu - Pm/t, 



dt 

dpR 
dt 



—tt = -4Hp R + pm/t, 



where H is the time-dependent Hubble parameter: 
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with p being the total energy density. (We assume a 
flat universe throughout). If the decaying nonrelativistic 
component dominates the expansion, it is easy to show 
that the energy density in the nonrelativistic component 
is rapidly converted into relativistic energy density when 
t ~ r (see, e.g., Ref. |7|); indeed, this assumption has 
become part of the standard cosmological lore. 

But now consider what happens for a vacuum- 
dominated expansion, such as the universe is experienc- 
ing at present. We take the ratio of radiation density to 
matter density to be given by r: 



r = Pr/p 
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Then equations (JTJ) — ((2|) can be combined to yield an 
equation for r: 



dr 
~dt 



--H 



(5) 



In the matter-dominated era, H decreases with time, 
and as long as H < 1/r, we see that r — > oo, as ex- 
pected. However, as the universe evolves from a matter- 
dominated state to a vacuum-energy dominated state, 
the value of H asymptotically approaches a constant 
value, Ha, given by 



H A = 



8nGpj* 
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where pa is the (constant) vacuum energy density. Define 
a 'time of no return', tA = ^/H\. Substituting H = 
Ha = l/t\ into equation (0, we can solve this equation 
analytically to yield 
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(7) 



where r is normalized to be zero at t = 0, i.e. the initial 
radiation density is negligible in the matter-dominated 
era. 

In terms of the present-day Hubble parameter Hq, and 
the fraction of the critical density in vacuum energy, ^a, 
we have simply 
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9.8 x 10 9 yr h^fl 
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where h is the value of Ho in units of 100 km sec 1 
Mpc" 1 . 

The asymptotic value of r is sensitive to the relative 
values of r and t\. If t < tA, then as t — > oo, equation ([7]) 
gives r — > oo, which is the same as the conventional re- 
sult for decays during the matter-dominated phase, since 
at such times the vacuum energy will indeed be cosmo- 
logically irrelevant. But if r > <a, the right-hand side of 
equation ([7]) goes to — 1 as t — > oo, and r asymptotically 
approaches a constant, given by 



r(t — » oo) = 



tA 



tA 



(9) 



Thus, the ratio of radiation to matter approaches a con- 
stant; the matter never disappears relative to the radi- 
ation (although of course both densities go to zero as 
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FIG. 1: The ratio of decay-produced relativistic energy den- 
sity, pr, to decaying nonrelativistic dark matter, pm, as a 
function of the time measured in units of the decaying parti- 
cle lifetime r. From top to bottom, the curves correspond to 
T/t A = 0.5, 1,2,5. 



t — > oo). Further, as long as r > 2£a 7 we have r < 1 
asymptotically; in this case the decay-produced radia- 
tion never dominates the decaying matter! Finally, there 
is an intermediate case, r = t\, for which the solution 
given in equation (0 does not apply. In the case, the 
solution of the evolution equation is simply r = tjr. In 
this case, r becomes arbitrarily large as t — > oo, but the 
increase is linear in time, rather than exponential as is 
the case for shorter particle lifetimes. (See also [l2| which 
examines a different model for energy exchange). 

While this result is perhaps non-intuitive, there is 
a straightforward physical explanation. Radiation red- 
shifts as one extra power of the scale factor relative to 
matter. For a vacuum energy-dominated universe the 
scale factor is itself an exponential of the time. Thus, 
in the absence of decay the energy density of matter in 
the universe will increase relative to radiation by an ex- 
ponential function of time. If the matter density itself 
decreases exponentially, the two factors cancel, leading 
to a constant final ratio of matter density to radiation 
density. 

Our analytic solutions give the correct asymptotic be- 
havior, but equation ([5]) must be integrated numerically 
to study the physically realistic case of a universe con- 
taining both decaying dark matter and vacuum energy. 
We have performed this numerical integration for several 
illustrative values of r/t\, (note that Pm/pa at any time 



is a fixed function of t/t\)', the results are displayed in 
Fig. 1. This figure clearly shows the asymptotic behavior 
discussed above. 

Now consider our present-day universe. Taking reason- 
able values of h = 0.7 and ft a = 0.7 in equation ©, we 
obtain a value of t\ — 1.7 x 10 10 yr. Current constraints 
on the decay of the dark matter yield r > t\. For ex- 
ample, Ichiki et al. [ll|, using the WMAP data, derive 
a 95% confidence limit of r > 5.2 x 10 10 yr, which cor- 
responds to r/t\ > 3. Thus, we conclude that we have 
already entered the era at which it is no longer possible 
for relativistic energy density from decaying dark matter 
to ever dominate the dark matter density itself. 

There are other mechanisms to convert the dark matter 
into relativistic energy density. For example dark mat- 
ter within the dense cores of halos can annihilate into 
relativistic particles (see, e.g., Ref. [l3l [l6f ). As Krauss 
and Starkman have recently demonstrated [l4| , such pro- 
cesses will have dramatic consequences for the future of 
large scale structure. 

In general, as described in at late times the an- 
nihilation rate for dark matter in halos with density n, 
given by Tann = n(av), will fall as t~ l . Nevertheless, 
even if it were to remain constant over time, this anni- 
hilation rate would still not be fast enough produce a 
radiation density that would overwhelm that in matter. 

In the case of falling annihilation rates one can in fact 
produce a stronger bound. If M — M/t, which will 
be the case for bound systems whose density changes 
only due to annihilations as the universe ages [l5[, then 
M(t) = Mo/t. In this case, as long as the expansion of 
the universe is such that R w t a where a > 1 then radia- 
tion will redshift faster than the mass decreases, and the 
radiation density will never overwhelm the matter den- 
sity. This will occur if the equation of state parameter 
w = P/p < — 1/3 for the dominant energy in the universe. 
However, for systems bound by gravity, where adiabatic 
expansion will accompany annihilations, supplanting an- 
nihilation in reducing core densities, then M = M/At 
[13 ]. In this case, M{t) = Mo/t 1 / 4 . For this case, no ex- 
pansion involving radiation or matter is sufficiently slow 
to allow the radiation produced by annihilations to over- 
whelm the matter density. 

In the unphysical but conservative case where annihi- 
lation occurs at a constant rate, dark matter annihila- 
tion would then map directly onto the decaying particle 
problem we have already solved. With a constant annihi- 
lation rate, each galaxy halo would act as a very massive 
decaying "particle", with lifetime r = 1 /Tann, and our 
previous analysis would then hold. Studies of galaxy halo 
profiles today (i.e. [HI, [l?]]) imply that the characteristic 
annihilation time, r, in galaxies must greatly exceed the 
present age of the universe, and hence also greatly ex- 
ceed t\ today. Moreover for annihilation into standard 
particles and for dark matter whose remnant abundance 
today was determined by freezeout at early times, the 
inverse annihilation rate will in fact be many orders of 
magnitude larger than the current age today. In either 
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case, our previous argument holds: the ratio of pa to Pm 
will remain small regardless of the particle annihilations. 

Our consideration of dark matter annihilation is par- 
ticularly relevant if ordinary matter is in fact unstable, 
due to proton decay. In this case, our arguments tell 
us that even if dark matter annihilates away, remnant 
baryons will still dominate compared to radiation. More- 
over, even once baryons start to decay significantly in 
the far, far future they will still dominate over the 
energy density of their relativistic decay products until 
these decay products themselves become non-relativistic, 
assuming they are massive. Thus, matter, even ordinary 
matter, will always continue to dominate over radiation 
for all times. 

We further point out in this regard that our analysis 
clearly applies in the case of two decaying matter species 
with different timescales, as for example might occur for 
decaying dark matter and decaying baryons. Since the 
latter are likely to have a far longer lifetime, our analy- 
sis is trivially extended by considering two species with 
different exponential timescales for decay. Consider two 
such species with lifetimes n and t%, with n << T2- 
Once t > n, while the ratio of energy densities in the 
first decaying particle and the radiation produced by it 
will become constant, nevertheless the total energy den- 
sity in both will quickly become insignificant relative to 
the energy density in the second, more stable, nonrel- 
ativistic species. Then when the second nonrelativistic 
species decays, our original analysis applies. 

In the extremely far future, individual galaxies, or clus- 
ters of galaxies will become essentially isolated 'island 
universes' |.20j. Nevertheless, on sufficiently large (per- 
haps even super-horizon) scales these will behave like a 
homogenous gas of particles amidst a background vac- 
uum energy sea from the perspective of this analysis, and 
if the matter within them decays or annihilates, as dis- 
cussed above, they can be treated as decaying particles 



and our analysis continues to hold. 

To be sure, it is true that in the far future of a vacuum- 
dominated universe, neither matter nor radiation will sig- 
nificantly affect the dynamics of the universe, and there- 
fore which dominates over the other is, from a practical 
perspective, not dynamically significant. Nevertheless, 
many scenarios in early universe cosmology are based on 
decaying systems which cause the universe to shift from 
matter to radiation domination. One's intuition suggests 
this will always be the case, and thus this result is inter- 
esting because it demonstrates that conventional wisdom 
about the past universe cannot be applied to the future 
universe, if it remains vacuum or dark energy dominated. 
In this regard, our result may also be relevant for deter- 
mining the initial conditions of a post- vacuum dominated 
universe, if the dark energy is not vacuum energy, but in- 
stead evolves. 

This latter comment, of course, reinforces the fact that 
our results depend crucially on the assumption that the 
dark energy is, indeed, a vacuum energy (cosmological 
constant) with pa = constant, and that it, itself, does 
not decrease with time, in which case the Hubble pa- 
rameter will once again decrease with time. Thus, in 
spite of the many negative facets of a future with a cos- 
mological constant [18| |. only persistent vacuum energy 
remarkably preserves this surprising eternal dominance 
of matter over radiation. 
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